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Abstract 

Rare b decay processes are analyzed in the multi-Higgs doublet model. Taking 
account of the constraint from the b — > sj process, the branching ratio and the 
forward-backward asymmetry of the final leptons for the b — > s£ + £~ process are 
calculated. It is shown that the branching ratio can be a few times larger than the 
standard model prediction and the asymmetry can be significantly different from 
that in the standard model. Combining these observable quantities it is possible 
to determine complex coupling constants associated with the charged Higgs mixing 
matrix. The CP violating charge asymmetry in the b — > s£ + £~ process and the 
branching ratio of the b — > s vv process are also calculated. 



In search for new physics beyond the Standard Model (SM), flavor changing 
neutral current (FCNC) processes can play an important role. Since the source of 
the flavor changing processes only lies in the Cabbibo-Kobayashi-Maskawa (CKM) 
matrix in the SM, it is possible to make precise predictions for various observable 
quantities in FCNC processes within the SM. An important example is the branching 
ratio of the 6 — > s 7 process, which was reported to be Br(6 ^37) = (2.32 ± 0.57 ± 
0.35) x 10~ 4 from the CLEO experiment [p]. Using the observed top quark mass, 
the measured branching ratio is consistent with the SM prediction. Therefore this 
process becomes a very strong constraint on models beyond the SM such as two 
Higgs doublet model and supersymmetric extensions of the SM || . In addition 
to the b — > s 7 process, the b — > s £ + £~ and b — > suu processes can be important 
constraints on new physics. Current experimental bounds for these processes are 
only by one order of magnitudes above the SM predictions [|, [|. 

In this letter we consider possible constraints on parameters in the multi-Higgs 
doublet model from these rare b decay processes. When the number of the Higgs 
doublets is more than three we can introduce complex phases in the charged Higgs 
mixing matrix. In this sense the multi-Higgs doublet model is a natural extension 
of the SM which involves a new source of the CP violation. Present phenomenolog- 
ical constraints on coupling constants in this model are considered, for example, in 
Ref. ||. Remarkably it is pointed out in Ref. that the strongest constraint on the 
imaginary part of coupling constant comes from the b — > S7 process, not from the 
CP violating quantities such as neutron electric dipole moment. Here we study the 
b — > s £ + £~ and b — > svv processes in the same model. We found that within the 
present experimental constraints including the b — > s 7 process the branching ratio 
can be a few times larger than the SM predictions and the forward-backward asym- 
metry of the final leptons for b — > s£ + £r process can significantly differ from the 
SM. Combining information from the branching ratio and the forward-backward 
asymmetry, it is possible to determine the complex coupling constants in future 
experiments. We also calculate the CP violating charge asymmetry between the 
b — > s£ + £~ and b — > s£ + £~ decays, which turns out to be a few percent in the 
allowed region. 

We consider the multi-Higgs doublet model with the following Yukawa cou- 
plings, 
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£ = QLVd d R H d + q L y u u R H u + £ L y t e R H £ + h.c, (1) 

where yd, y u , ye are 3x3 matrices. We have assumed that up-type quarks, down-type 
quarks and leptons have Yukawa couplings with only one Higgs doublet, H u , H d , Hi 
respectively. In this way we can avoid FCNC effects at the tree level ||. The 
couplings of physical charged Higgs bosons with fermions are then given by 

n-1 

C = (2V2G F ) 1/2 £ (Xi u L V M D d R Ht + Y, u R M v V d L H+ + Z { V L M E e R H? ) , (2) 
i=i 

where n is the number of Higgs doublets, Hf{i — 1 ~ n — 1) represent mass eigen 
states of charged Higgs bosons and V is the CKM matrix. New CP violating complex 
phases arise in the charged Higgs mixing matrix if there are three or more Higgs 
doublets. In such cases the coupling constants Xi, Yi, Zi are in general complex 
numbers. These are several relations among Xi, Yi, Zi from the requirement of 
unitarity of the mixing matrix . 

The charged Higgs interactions in Eq.@ can induce extra contribution to 
FCNC processes in this model. Here, we are interested in FCNC processes related 
to the b quark. Inclusive branching ratios of b —>■ 57, b — > s£ + £~ and b — > s vv are 
calculated through the following weak effective Hamiltonian at the bottom scale as 
described in [g 0; 

AC 10 

i=i 



The relevant operators at the bottom scale are 



7 = j^m b s L a^b R F^, (4) 



Os = -^- 2 m b s L T a a^b R G a ^ u , (5) 

O9 = ^s L rb L lj/, (6) 

010 = -^fSL^b L J ltil H, (7) 

011 = 1R f- 2n SL-fbL E %(l-7Vi, (8) 

167r 2 sm 6 W i= ^ )T 
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where e, g s are QED and strong coupling constants respectively, and 9w is the weak 
mixing angle. Using renormalization group equations of QCD the Wilson coefficients 
Ci's at the bottom scale are related to those at the weak scale. New physics effects 
enter through the Wilson coefficients at the weak scale. From the charged Higgs 
interactions these coefficients receive the following new contributions, 



1 n ~ 1 r /2 \ 



(9) 



+Xi Y* [ -F 3 (x thi ) + F 4 (x thi ) 



n-l 



°10 



-« E x tK m'F^) + X t Y*F 3 (x tht ) 
A i=i 



sin ^14/ 



C 



sin 2 # 



(10) 

(11) 

(12) 
(13) 



Here C H and are given by 

n-l i 



1=1 

n-l 



/2 

i=i v,:) 



(14) 



(15) 



where a^vv = m 2 / m wi x thi = m t/ m2 Hi- The functions Fi-Fq are defined as follows: 
Fi(x) = 12(:r 1 „ 1) 4 {x 3 - Qx 2 + 3x + 2 + 6xlnx), F 2 (x) = 12( ^ 1)4 (2a; 3 + 3a; 2 - 6a; + 1 - 



6a; In a;), -Fsfrr) 



2(x-l) 3 



4x + 3 + 2 In a;), FAx) 



2{x— i) 3 ("^ 1 2a; In a;), 
1 ' 11a; 3 + 



/•-,(.' ) = (7a; 3 - 36a; 2 + 45a; - 16 + (18a; - 12) lux), F 6 (a;) = 

18a; 2 — 9a; + 2 + 6a; 3 In a;). The Wilson coefficients at the weak scale are then given 
by 



C i (m w )=C? M {m w )+C?{m w ), 



(16) 
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where Cf M is contribution from the SM |], After taking account of the QCD 
corrections at the leading logarithmic order, Ci(m b ) can be expressed as, 

-7/i)+CV, (17) 

(18) 
(19) 
(20) 
(21) 

where r\ = a s (mw) / 'oc s (m b ) . CV ~ C 9 are constants which depend on the QCD 



C 7 (m b ) 


16 . s O 

= C 7 (m w )r]23 + C g (m ff )-(ij 


C 8 (m b ) 


= C 8 (m H /)^i + C 8 , 


C 9 (m b ) 


= C 9 (m w )+C g , 


C w (m b ) 


= C 10 {m w ), 


C n (m b ) 


= C xl {m w ), 



coupling constant. Detailed formulas are found in [JTOj, [II]]. Numerically these are 
given by C 7 = -0.17, C 8 = -0.077, C 9 = 1.9 for a s (m z ) = 0.12. 
The b — > s 7 branching ratio is then given by 



Br(6 — ► 57) = Br(6 -tcev) 



cb 



\C 7 (m b ) 



(22) 



7rg(m c /m 6 ) 

where the phase space factor g(z) is given by g(z) = 1 — 8z 2 + 8z 6 — z 8 — 24z A \nz. For 
the b — > s£ + £~ process, the differential branching ratio and the forward-backward 
asymmetry of leptons in the lepton-center-of-mass frame are given by 



dBr(b -> s£+t 
d s 



Br (6 — > ceu)- 



4tc 



v ts v t t 2 



V, 



cb 



g{m c /m b ) 



(23) 



(|C 9 + Y(s)\ 2 + \C w \ 2 )(l + 2s) + -\C 7 \\2 + s) 2 + 12ReC 7 *(C 9 + Y(s)) 

s 



A FB (s) 



(24) 



-3ReC 1 * [(a 9 + y(s))s + 2C 7 



~ (|C 9 + F(s)|2 + |C 10 | 2 )(1 + 2s) + l\C 7 \ 2 (2 + s) 2 + 12ReC 7 *(C 9 + F(s)) ' 
where s = s/m 2 , = (p++p~) 2 /m 2 . andp + Qo_) is the four momentum of £ + (£~). Y(s) 



represents the contribution from the charm quark loop at the m b scale. See Ref . JTT 
for details. Here we neglect the J / ip and ip' resonance contributions. The branching 
ratio for the b — > s vv is given by 



Br(6 — > s ViVi) = 3 Bi(b — > c e v) 



a 



47r 2 sin^ 



VtbV* 



V, 



cb 



g(m c /m b ) 



(25) 
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Using the above formulas, it is now straight-forward to evaluate the branching 
ratios and the asymmetry numerically in the multi-Higgs doublet model. In the 
followings, we assume that only one of the physical charged Higgs boson is light 
and neglect the effects of other physical charged Higgs bosons. Then we keep only 
one term in the summation in the expressions for Cf's. Dropping the index i for 
the lightest charged Higgs bosons, the relevant parameters are \Y\ 2 , XY* and the 
mass of the charged Higgs boson. It should be noted that the coefficients CV and C$ 
depend on both \Y\ 2 and XY* whereas C 9 and C w only contain \Y\ 2 . 

Since the b —>■ s'j branching ratio is already observed and is consistent with 
the SM prediction, we use Br(6 — > 57) to solve ImXy* in terms of \Y\ 2 and KeXY* 
for each value of the charged Higgs mass. Then other observable quantities for 
b — > s £ + £~ and b — > s vv processes can be calculated as functions of two parameters. 
We should also take into account constraints from other processes. Besides the 
b — > 57 process, the B°-B° mixing and the Z — > bb process give the most strong 
constraints on the possible value of \Y\ 2 |J. The contribution to the mixing 
from the charged Higgs boson is expressed as 

Mg = ^y w r, B {V u V; b ) 2 ^B B f 2 B m B x t [yMyt)\y\ A (26) 

+x t (2I 2 (x t ,x H ) -8h(x t ,x H )\Y\ 2 ) 

where yt = m 2 /m 2 H , x B = ™ 2 H/ m wi an d the functions /1-/3 are defined as follows. 

r / \ 1+ x 2x\nx /ri _. 

/iW = (T^w + jr^' (27) 

x y 2 \n(y) x(-x - xy + 2y) ln(x) 

2[X,V) (x-y)(x-l) (y-l)(x-y) 2 {I - x) 2 {x - y) 2 ' 1 ' 



u x y \ = 1 + y ln y + (-^ 2 + y) M x ) ^ 

(x — y)(x — 1) (y — l)(x — y) 2 (l — x) 2 (x — y) 2 



Here we have retained only the relevant terms which are proportional to \Y\ and 
\Y\ 2 . Since this quantity depends on the CKM matrix element V t d which has not 
been known well, the constraint from B^-B^ mixing is not very strong. Using the 
constraint on Vtd from the charmless b decay (0.005 < \Vtd\ < 0.012) and taking 
account of uncertainties from f B and B B (f B y/B B = 200 ± 40MeV) we can deduce 
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\Y\ < 1.5(2.1) for m H = 100(300) GeVQ. Here and in the followings we fix the top 
quark mass as 175 GeV. For the Z — > bb process we calculated the charged Higgs 



contribution to R b = T(Z — > bb)/T(Z — » hadrons) following Ref. ||12||. We define 
the deviation from the SM contribution (i?f M = 0.2158) by 5R b = R b - R$ M . The 
charged Higgs contribution to 5R b is shown in Fig.|l] as a function of \Y\ 2 . In this 
figure we also show the 3 a lower bound of 5R b from the world average value R, 



b 



0.2178 ± 0.0011 @. From Fig.| we can derive \Y\ < 0.9(1.3) for m H = 100(300) 



GeV. Note that if we use the value R b = 0.2159±0.0009(stat)±0.0011(s|/st) reported 



by ALEPH []14| the lower bound of 5R b shifts to —0.0041, which corresponds to the 
upper bound \Y\ < 1.6(2.3) for m# = 100(300) GeV. Since experimental situation 
is not conclusive, to be conservative, we present the result in a rather wider range 
of \Y\. 

Fig.Q(a) shows that the ratio of Br (6 — > s /i + /U~) branching ratio for the multi- 
Higgs doublet model normalized by the SM branching ratio in the space of \Y\ 2 
and Ke(XY*) for uih = 100 GeV. For numerical calculations we use m b = 4.7 GeV 
and m c = 1.5 GeV. To avoid the large effect of the J/ip resonance we integrate the 
differential branching ratio for the kinematical range 4m^ < s < {rrij/^p — 5) 2 where 
is muon mass, rrij/^p is J/ip mass and 5 = 100 MeV. The interference effect 
between the J/if> resonance and the short distance contribution in this kinematical 



region is still sizable (< 20%) for the SM case [fLlj| . However as we can see in Fig.g the 
interference effect is expected to be smaller than the charged Higgs contribution to 
the short distance part. In order to solve lm(XY*) in terms of \Y\ 2 and Ke(XY*), we 
assume Br(6 — ► sj) = 2.8 x 10 -4 . Since there is a sizable experimental error on this 
quantity we have just used the SM prediction for the illustration. We expect that the 
experimental error as well as the theoretical ambiguity will be reduced in future when 
actual analysis on b — > s£ + £~ is done. Fig.^(b) shows the asymmetry in the same 
kinematical range. From these two figures we can see that the branching ratio can be 
a few times lager than that of the SM, which is 3.8 x 10 -6 in this kinematical region, 
and the asymmetry can be as large as 20 % compared to about 5 % in the SM. It is 

* The CP violating parameter in the K°-K° mixing, ex, also receives a similar contribution 
from the charged Higgs loop as the B°-B mixing. In this case, however, the relevant CKM matrix 
element is different and the constraint on |V| from ck only is not strong. If we combine the B°-B 
mixing and the 6k constraints we can exclude a slightly larger parameter space but the allowed 
region of \Y\ is numerically almost the same as above. 
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also interesting to see that the branching ratio is sensitive to the value of \Y | 2 , on the 
other hand the asymmetry gives an independent constraint on the parameter space. 
Therefore the values of \Y\ 2 and Ke(XY*) are determined if the branching ratio and 
the asymmetry are measured with reasonable accuracy. Fig||(c) shows the b — > s W 
branching ratio normalized by the SM prediction (Br(6 — > svT>)sm — 4.3 x 10~ 5 ) 
in the same parameter space. In this case the branching ratio only depends on 
the parameter \Y\ 2 and is enhanced by a similar factor as the branching ratio of 
b^s£+£~. 

If lm(XY*) ^ 0, there is possibility to observe the CP violating charge asym- 
metry i.e. the difference of b — > s £ + £~ and b — > s£ + £~ . Since this quantity is induced 
with help of the phase in Y(s) in Eq. (|23|) , the asymmetry appears only above the 
cc threshold of the lepton invariant mass. The charge asymmetry is defined as 

[■ml , ( dBr{b^s l+t~) _ dBrjb^s 1+ i~) \ 

J T ^{ » » 

« j / dBr(b^s£+£-) , rfBr(b->sl+l-) \ 

km^+S) 2 0,8 I ds + ds J 

where m^i is the ip' mass and we take 5 = 100 MeV. In Fig.||(d) this quantity is shown 
for rriH = 100 GeV. Since this asymmetry is at most a few percent, we need large 
statistics to measure it. Note that there is a twofold ambiguity in sign of lm(XY*). 
FigfJ(a)-(c) do not depend on the sign of lm(XY*) while Aqp (Fig.[|(d)) change its 
sign according to the sign of \m.[XY*\ Fig.|2|(d) shows the case of lm(XY*) > 0. 

In Fig.^(a)-(d) we show the same quantities for m# = 300 GeV. As in the case 
of mil — 100 GeV we can see similar enhancements on the branching ratios and the 
forward-backward asymmetry in the allowed parameter space. 

In summary we have investigated possible constraints on the charged Higgs 
coupling parameters from the b — > sj, b — > s£ + £~ and b — > svv processes. Within 
the present experimental constraints including 6—^57 process, we have shown that 
the branching ratio of the b s i + £~ and b —>■ s vv processes can be enhanced 
a few times compared with the SM. The lepton forward-backward asymmetry in 
the b —>■ s£ + l~ process can also be enhanced by a factor of several. Since these 
quantities depend on the model parameters in different ways, we can obtain useful 
information on the model once these processes are observed experimentally. 

This work is supported by the Grant-in-aid for Scientific Research from the 
Ministry of Education, Science and Culture of Japan. 
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FIG. 1: The charged Higgs contribution to 5Rt> for m H = 100 GeV(solid line) and 
ran = 300 GeV(dashed line). The dotted line represents the 3 o lower bound from 
the world average value of Rf,. 



FIG. 2: (a). The ratio of the Br(6 — > s£ + £~) in the range 4m^ < s < {mj/^ — 5) 2 
for the multi-Higgs doublet model (MHDM) and the standard model(SM) in the 
case of ma = 100 GeV. (b).The forward-backward asymmetry {Aps) in the range 

Ami - s - ( m J/Tp ~ 5 ) 2 for the MHDM and tne SM in tn e case of m H = 100 
GeV. (c).The ratio of the Br(6 svv) for the MHDM and the SM in the case of 
mH = 100 GeV. (d).The charge asymmetry (Acp) in the range (m^i +5) 2 < s < ml 
for m H = 100 GeV. We take 5 = 100 MeV in (a), (b) and (d). 



FIG. 3: The same as Fig. | for m H = 300 GeV. 
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